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ABSTRACT 



Context. The mechanism by which supergiant (sg)B[e] stars support cool, dense dusty discs/tori and their physical relationship with 
other evolved, massive stars such as luminous blue variables is uncertain. 

Aims. In order to investigate both issues we have analysed the long term behaviour of the canonical sgB[e] star LHA 115-S 18. 
Methods. We employed the OGLE II-IV lightcurve to search for (a-)periodic variability and supplemented these data with new and 
, , ' historic spectroscopy. 

', Results. In contrast to historical expectations for sgB[e] stars, S18 is both photometrically and spectroscopically highly variable. 

The lightcurve is characterised by rapid aperiodic 'flaring' throughout the 16 years of observations. Changes in the high excitation 
emission line component of the spectrum imply evolution in the stellar temperature - as expected for luminous blue variables - 
although somewhat surprisingly, spectroscopic and photometric variability appears not to be correlated. Characterised by emission in 
low excitation metallic species, the cool circumstellar torus appears largely unaffected by this behaviour. Finally, in conjunction with 
intense, highly variable He n emission, X-ray emission implies the presence of an unseen binary companion. 

Conclusions. S18 provides observational support for the putative physical association of (a subset of) sgB[e] stars and luminous 
,<_» ' blue variables. Given the nature of the circumstellar environment of S18 and that luminous blue variables have been suggested as SN 

progenitors, it is tempting to draw a parallel to the progenitors of SN1987A and the unusual transient SN2009ip. Moreover the likely 
binary nature of S18 strengthens the possibility that the dusty discs/tori that characterise sgB[e] stars are the result of binary -driven 
mass-loss; consequently such stars may provide a window on the short lived phase of mass-transfer in massive compact binaries. 

Key words. stars:evolution - stars:early type - stars:binary - star:individual:LHA 115-S 18 - star:emission-line 

On 

in 

t^J- | 1. Introduction I1995I Maeder & Mevnet l200T1 l. there is no guarantee that the 

d^i . B[e] phenomenon is attributable to a single physical mechanism, 

• , Supergiant (sg)B[e] stars are thought to represent a short lived since stars at the extremes of this d i stri bution will follow very 

■ phase in the evolution of massive stars from the main sequence different post . main seqU ence evolutionary paths. Indeed, the in- 

^ ; through to H-depleted Wolf-Rayets (WR). Observationally, they dusion of sgB[e] stars m such schemes is highly uncer tain at 

are characterised by an hybrid spectrum consisting of broad p rese nt (Aret et al 1 201 2t 
emission lines in high excitation species and narrow, low excita- 
tion metallic lines and a near-mid IR continuum excess. Zickgraf 

et al. d 19851 1 proposed a two component circumstellar environ- Despite their co-location in the Hertzsprung-Russell (HR) 
ment to reconcile these findings, comprising a hot, high velocity diagram (cf. Clark et al. 2005b. Aret et al. 2012), it is cur- 
I polar wind and a low velocity equatorial disc or torus. The latter rently unknown whether sgB[e] stars and luminous blue vari- 
is expected to be either slowly expanding or in quasi- Keplerian ables (LBVs) may be physically identified with one another. As 
rotation (Zickgraf et al 120031 Kraus et al. 1201 01 Aret et al. 120 121 their name suggests, LBVs are massive evolved stars that are 
and appears to be cool an d dense en ough to permit both neutral subject to pulsational instabilities that drive dramatic photomet- 
ga s (Zick graf et al. 1 19891 Kraus et al.|2007) and dust (Kastner et ric and spectroscopic variability (e.g. Humphreys & Davidson 
al. 120061) to exist. 119941 Sect. 4.1) leading to changes in stellar radius, tempera- 
While this model is successful in explaining the observa- ture and in some cases, mass loss rate and bolometric luminos- 
tions, there is currently no consensus as to either the detailed ge- ity (e.g. Groh et al. 2009 Clark et al. 2009b and refs. therein), 
ometry of the circumstellar environment or the physical mecha- Conversely, sgB[e] stars have been considered to demonstrate 
nism by which it could form. With bolometric luminosities rang- little photometric variability (e.g. Zickgraf et al. 1986), with 
ing from ~ 10 4 - 1O 6 L implying intial masses ranging from only a handful of known counter-examples (Sect. 4.1). Similarly, 
1OM to > 6OM (Zickgraf et al. 119861 Gummersbach et al. spectroscopic variability consistent with the canonical LBV 

'outbursts' amongst sgB[e] stars has yet to be observed, although 

* This work is partly based on observations collected at the European both classes of objects may show c onside rable overlap in spec- 
Southern Observatory (programme ID ESO 088.D-0352(A). tral morphologies (e.g. Morris et al. [1996). 
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Fig. 1. Top panel: plot o f the 1997-2012 OGLE-II, -III and -IV lightcurve data for LHA 1 15-S 18 with the timings of optical blue- 
(Massey & Duffy 2001, Torres et al. 2012 and this work) and far-red (Aret et al. 120121 and this work) spectroscopic and X-ray 
(Antoniou et al. 2009 Bartlett et al. in prep.) observations indicated. The nearest photometric data points to each observation are 
given by the filled red symbols. Note that the spectrum of Massey & Duffy and the first epoch of Torres et al. were obtained quasi- 
simultaneously and the date indicated for the first X-ray observation is the midpoint of the period over which observations were 
made. Photometric errors are smaller than the symbol size plotted. Lower panels: expanded presentation of two sub-sections of the 
lightcurve (indicated by dotted lines in the upper panel). 



In order to address these issues we present new observations 
of the sgB[e] LHA 115-S 18 (henceforth S18) which has been 
advanced as an LBV by a number of authors (e.g. Morris et al. 
[19961 van Genderen & Sterken|2002}. 

Since its discovery by Henize (1956), numerous observa- 
tions have been made of SI 8. The optical spectrum of the star 
is characterised by strong emission in the Balmer series and 
prominent He i lines (e.g. 6678A and 7065 A) as well as a wealth 
of forbidden and permitted low excitation metallic lines (Fe n, 
[Fen], Tin, Cm and Oi). Following the behaviour observed in 
other sgB[e] stars, the high excitation lines are broader than the 
metallic lines. Moreover, they exhibit pronounced variability in 
profile and strength - the Balmer lines evolve from (asymmetric) 
single peaked to P Cygni profiles with absorption extending to 
750kms~ 1 and vary in strength by a factor of ~ 3 (Zickgraf et al. 
1989 Nota et al. 119961 see Table A.l for a historical summary 
of optical spectroscopic studies). Recently, high resolution ob- 
servations by Torres et al. (12012) revealed that the metallic lines 



are also variable, while Nota et al. (1996) reported changes in the 
radial velocities of these lines. However the most dramatic vari- 
ablity is seen in the He n 4686A line, which has been observed 
to transition from absence to being comparable in strength to H/3 
(Table A.l), with such changes occuring over short timescales 
(of the order of months; Shore et al. 11987) . 

Contrary to the conclusions of Shore et al. (1987) and 
Zickgraf et al. (fT989l . van Genderen & Sterken (120021 found 
S18 to be photometrically variable at optical wavelengths be- 
tween 1987-1991. Modulation occured on timescales from days 
to years, with peak-to-peak amplitudes ranging from ~ 0.1 - 
1.0 mag (see Sect. 2.1.2. for a full description). 

At other wavelengths, five International Ultraviolet Explorer 
UV observations between 1981 July to 1983 March were re- 
ported in Shore et al. (I1987I I. No photometric variability above 
~0.25mag. was found but a number of high excitation wind fea- 
tures (Hen, Civ, Nv) were significantly variable. Two epochs 
of I band spectroscopy made by Shore et al. (1987 obtained 
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Fig. 2. Upper panel: un-detrended Lomb-Scargle periodogram 
for the I-band photometric data, with the position of 1/year indi- 
cated by the dashed line. No periodic signal is returned. Lower 
panel: lightcurve with a sine curve corresponding to the ~440 d 
period implied by the largest peak in the periodogram overplot- 
ted - the periodicity is not stable throughout the 16 year duration 
of observations. 

on 1983 October 10) and Aret et al. (120T21 obtained on 2005 
December 10-11) share the same bulk morphology; being dom- 
inated by the Paschen series, Can 8491, 8542 and 8662 A 
and O i 8446A. Conversely, K band spectroscopy presented by 
McGr egor e t al. (119891 obtained between 1987-89) and Morris 
et al. ( 11996t obtained in 1995 November) reveal evolution over 
this period, with the development of He i 2.058yum and 2. 1 12//m 
and CO bandhead emission between the two epochs. At longer 
wavelengths, mid-IR spectroscopy and photometry provided by 
Kastner et al. (120101) and Bonanos et al. ( 120101) show that S18 
supports a substantial IR-excess. 

Finally a single epoch of Chandra X-ray observations made 
between 2002 May-October (source 6_20 of Antoniou et al. 
2009) provides a weak detection of SI 8. While the count rate 
was too low to model the spectrum, they estimate a flux of 
~ 3 x 10 33 erg s _I from this observation. Bartlett et al. (in prep.) 
present the results of an X-ray survey of sgB[e] stars within the 
Magellanic Clouds and provide a confirmatory XMM-Newton 
detection from an observation on 2003 December 18. Further 
details on data reduction and analysis are contained therein, but 
utilising the count rate of the EPIC-pn detector and Webpimm^] 

1 Assuming f? H = 6x 10 20 cirr 2 (corresponding to the column density 
in the direction of the SMC; Dickey & Lockman 19901, T = 1.7 (similar 



it was possible to obtain a preliminary luminosity estimate of 
~ 3.9 ± 0.7 X 10 33 erg s -1 . 

The remainder of the paper is structured as follows: we 
present our new photometric and spectroscopic observations, 
data reduction and analysis in Sects. 2 and 3. In Sects. 4 and 
5 we discuss the nature of S 18 in light of these results and possi- 
ble physical associations between sgB[e] stars, massive binaries 
and LB Vs. Finally we summarise our findings in Sect. 6. 



2. Data Reduction 

2.1. Photometry 

We utilized the Optical Gravitational Lensing Experiment data 
(OGLE-II, -III and -IV, Udalski et al. [19971 Udalski [20031) to 
investigate the long-term behaviour of S18 (OGLE-II designa- 
tion: SMC_SC6 311 169). OGLE has been regularly monitoring 
this object since 1997 with the 1.3-m Warsaw telescope at Las 
Campanas Observatory, Chile, equipped with three generations 
of CCD camera: a single 2kx 2k chip operated in driftscan mode 
(OGLE-II), an eight chip 64 Mpixel mosaic (OGLE-III) and fi- 
nally a 32-chip 256 Mpixel mosaic (OGLE-IV). Observations 
were collected in the standard V- and I-bands with the vast ma- 
jority of images taken through the I-band filter. 

OGLE I-band observations of S18 required an individually 
tailored approach to data reduction. During the period of OGLE- 
II observations the stellar brightness was well within the pho- 
tometic range of the OGLE-II camera and measurements were 
derived with the standard OGLE DIA photometric reduction, 
based on the image difference technique (Szymanski 2005). 
However, during the OGLE-III and OGLE-IV phases the bright- 
ness of S 18 on the survey frames turned out to be close to the 
saturation level of the CCDs. In particular this happened during 
observations with very good seeing conditions and as a result a 
few central pixels of the star saturated. 

Therefore in the case of OGLE-III observations we decided 
to employ point spread function (PSF) photometry obtained with 
the DoPhot photometry program (Schechter et al. 1993 also pro- 
vided by the OGLE-III photometry pipeline). The PSF approach 
is well suited to the treatment of saturated pixels at the cost of 
somewhat larger photometric errors. In the case of OGLE-IV 
images we manually derived S18 photometry, limiting ourselves 
to the images taken during seeing conditions where the object 
was not saturated. The photometric technique was identical to 
the OGLE standard photometric pipeline (Udalski 12003) . 

V-band images collected at all OGLE phases were reduced 
with the standard OGLE DIA technique. S18 was far from the 
saturation level in all OGLE V-band frames. 

All OGLE datasets were then combined after applying sec- 
ond order corrections for slightly different filter passbands. The 
complete dataset comprises 1365 I-band, 133 V-band and 100 
(V-I) colour index measurements over the 16 year period from 
1997-2012. Given the rapid variability present in the lightcurve 
we required both V- and I-band measurements on the same night 
to derive the (V-I) colour index. Errors on the OGLE DIA pho- 
tometry (I-band: OGLE-II and OGLE-IV, V-band) are well be- 
low 10 mmag, while those for PSF photometry (I-band: OGLE- 
III) below 20 mmag. Given the much better sampled I-band 
lightcurve we base our discussion on this, utilising the V-band 
data to provide additional colour information. 



to that adopted by the 2XMM survey; Watson et al.2009 ) and a distance 
of ~ 61 ± 3kpc (Hilditch et al. l2005i 
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Fig. 3. Plots of V- and I-band OGLE lightcurves (with the V-band displaced by 1.2 mag. to enable easy comparison) and (V-I) 
colour index of time. Timings of optical spectroscopic observations indicated by arrows and (V-I) data within 9 d of an observation 
indicated by the red symbols (with the first four data points within 3 d). 
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Fig. 4. Colour magnitude plot for S 1 8 . 
2.2. Spectroscopy 

The spectroscopic data presented in this work come from a num- 
ber of sources. Firstly, we reproduce the full spectrum from 
Graus et al. (120121 1. This was obtained on 2010 November 24 
with the MIKE echelle spectrograph on the Magellan Clay tele- 
scope and has a resolution of ~ 28, 000 over the 3250A-5050A 
wavelength range, with details on reduction provided in Graus et 
al. (120121 1. 



Observations taken on 201 1 July 20th and 21st and 2012 July 
7th were made with the fibre-fed dual-beam AAOmega spectro- 
graph on the 3.9 m Anglo- Australian Telescope (AAT) at the 
Australian Astronomical Observatory. The instrument was op- 
erated with the Two Degree Field ("2dF") multi-object system 
as front-end. Light is collected through an optical fibre with a 
projected diameter of 2'.'1 on the sky and fed into the two arms 
via a dichroic beam-splitter with crossover at 5700A. Each arm 
of the AAOmega system is equipped with a 2kx4k E2V CCD 
detector and an AA02 CCD controller. The blue arm CCD is 
thinned for improved blue response. The red arm CCD is a low- 
fringing type. On both occasions, the red arm was equipped with 
the 1700D grating, blazed at 10000A. This grating provides a 
resolving power R = 10000 over slightly more than 400A. In 
the 201 1 run, the central wavelength was set at 8600A. In 2012, 
it was centred on 8690A. The exact wavelength range observed 
depends on the position of the target in the 2dF field. 

During the 2011 run, the blue arm was equipped with the 
1500V grating, which provides R = 3700 over ~ 750A. The 
central wavelength chosen was 4400A. During the 2012 run, we 
used grating 580V, giving R = 1300 over ~ 2100A. The central 
wavelength was set at 4800A. 

Data reduction was performed using the standard automatic 
reduction pipeline 2dfdr as provided by the AAT at the time. 
Wavelength calibration was achieved with the observation of arc 
lamp spectra immediately before each target field. The lamps 
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provide a complex spectrum of He+CuAr+FeAr+ThAr+CuNe. 
The arc line lists were revised and only those lines actually de- 
tected were given as input for 2df dr. This resulted in very good 
wavelength solutions, with rms always < 0.1 pixels. 

Sky subtraction was carried out by means of a mean sky 
spectrum, obtained by averaging the spectra of 30 fibres located 
at known blank locations. With this procedure, any existing neb- 
ular emission is not removed. 

The optical data from 2011 Decemeber 8 were taken with 
the ESO Faint Object Spectrograph (EFOSC2) mounted at 
the Nasmyth B focus of the 3.6m New Technology Telescope 
(NTT), La Silla. The EFOSC2 detector (CCD#40) is a 
Loral/Lesser, Thinned, AR coated, UV flooded, MPP chip with 
2048x2048 pixels corresponding to 4. 1'x4. l'on the sky. The in- 
strument was in longslit mode with a slit width of 1.5"using 
Grism 14 - resulting in a wavelength range of AA3Q95 - 5085 A - 
and a grating of 600 lines mm" 1 and a dispersion of 1 A pixel" 1 
yielding a wavelength range of AA60Q0 - 7150 A. The result- 
ing spectra have a spectral resolution of ~ 10 A. The data 
were reduced using the standard packages available in the Image 
Reduction and Analysis Facility (IRAF). Wavelength calibra- 
tion was implemented using comparison spectra of Helium and 
Argon lamps taken through out the observing run with the same 
instrument configuration. 



3. Data presentation and interpretation 

3.1. Photometry 

Photometrically, S18 is clearly highly variable on multiple 
timescales (Fig. 1), mirroring the findings of van Genderen & 
Sterken (2002) and raising the possibility that such behaviour 
has continued uninterrupted throughout the ~ 22 year period fol- 
lowing their 1987-1991 observations. Given this, we followed 
the technique outlined in Bird et al. (2012) to search for peri- 
odic modulations in the I-band lightcurve. Briefly, we used the 
fast implementation of the Lomb-Scargle periodogram (Press 
& Rybicki 119891 ) to search for periodicities between 2 d (the 
Nyquist frequency for daily sampled lightcurves) and 10 3 d 
(above which the finite length of the lightcurve compromised 
sensitivity). Given the difficulties in detecting high-frequency 
signals in the presence of low-frequency modulations we also at- 
tempted to de-trend the data. Moving in one day steps, a rolling- 
mean was calculated over a ~ 50 - 300 d time window to de- 
termine a detrending model which was then subtracted from the 
light curve. 

No periodicity was identified (Fig 2). Despite the appearance 
of various subsets of the lightcurve (e.g. JDs~ 2455200 - 6100) 
no quasi-period present is stable over the 16 year duration of 
the observations and hence we are forced to procede via a more 
qualitative analysis, van Genderen & Sterken (2002) identified 
3 types of variability: (i) short-term variability on the timescale 
of days and ~ 0.1 — 0.2 mag amplitude, (ii) ~ 150 d modulation 
with amplitude ~ 0.25 - 1 .0 mag. and (iii) a long term (~ 2 year) 
trend with amplitude ~ 0.7 mag upon which the other behaviour 
is superimposed. We find no evidence for the latter behaviour 
in our data, seeing no long-term trends in brightness, but can 
recognise the other variability described. 

Specifically, we find S18 to vary between I~ 12.1 - 13.3 mag 
over the course of our observations but cannot identify any pro- 
longed 'quiescent' maxima or minima in this period. Rapid, 
low amplitude (~ 0.2 mag) modulation on timescales of days- 
weeks is present - e.g. the periods between JD~2453500-800, 



2453900-4150 and 2455250-5600 (Fig. 1). This may be super- 
imposed on the larger amplitude fluctuations that correspond to 
the second variability type of van Genderen & Sterken (12002) - 
e.g. JD-2452 100-300 and 2454700-800. The latter events span 
a range of timescales and amplitudes and appear exemplified by 
the behaviour of S18 between JD-2454000-4800 (Fig. 1). The 
initial fade of ~ 1 mag amplitude - upon which is superimposed 
a short-lived ~ 0.2 mag flare - occurs over a period of 400 d and 
is followed by a ~ 0.85 mag brightening event of only ~ 50 d 
duration. Note that we do not claim that the two types of photo- 
metric variability are the result of physically distinct processes 
within the system and merely employ them as a qualitative de- 
scription of the lightcurve; indeed it is entirely possible they may 
both be manifestations of the same underlying instability. 

Interpretation of colour variability is likewise difficult and 
is complicated by the comparatively sparse and uneven sam- 
pling of the V-band data, van Genderen & Sterken (2002]) re- 
port that at the photometric maxima corresponding to ~150 d 
variations the system appears bluer and we find such a broad 
description to be applicable to our data; the magnitude of vari- 
ability appears larger in the V- rather than the I-band (Fig. 3). 
However this behaviour appears not to be simply proportional to 
the brightness of the system as is clear from the colour mag- 
nitude plot (Fig. 4); such a relationship is only apparent for 
V< 13.6 - 13.8 mag and below this variability occurs at con- 
stant colour. Indeed, the 'blue-ing' of the system appears to be 
associated with episodes of brightening, with the subsequent 
decay branch of the lightcurve occuring at redder colours ((V- 
I)~1.2) - e.g. the 'flaring' episodes between JD-2450600-800 
and 2455700-6000. Unfortunately, better sampled data will be 
required to fully explore these relationships. 

3.2. Spectroscopy 

Globally, the four new blue-end spectra obtained between 2010 
November and 2012 July fit seemlessly into the long-term pat- 
tern observed for S18 (Table A.l and refs. therein), although 
the high resolution and signal to noise (S/N) ratio of three of 
the spectra enable us to identify qualitatively new behaviour. 
All four epochs are dominated by Balmer series emission lines, 
which demonstrated P Cygni profiles in all transitions in 2010 
November and Hy and higher lines in 2011 December (Figs. 
5-7fl Strong Hen 4686A emission was present in two spec- 
tra (2011 and 2012 July), weakly in emission in a third (2011 
December) and completely absent in our highest S/N and resolu- 
tion spectrum (2010 November). We highlight that the presence 
of strong He n 4686A emission appears anti-correlated with P- 
Cygni profiles in the Balmer series in these data, with the 201 1 
December spectrum apparently representing a transitional case. 
Similar behaviour is visible in the spectra from 2001 and 2005- 
2008 (Torres et al. 120121 ). and also appears present in earlier ob- 
servations (Table A.l). 

The appearance of Hen also heralds an increase in the 
strength of the higher Balmer series lines (Fig. 6) and we are also 
able to identify a corresponding increase in strength of weaker 
He i and He n lines which had hitherto been lost in the forest of 
low excitation metallic lines. Of particular interest are the ap- 
pearance of N m 4634A and 4640A emission lines that also ac- 
company Hen 4686A (Fig. 5). Conversely, while the wealth of 

2 Velocities of ~500±10kms~' for the blue edge of the absorption 
component in H/? and Hy were found for the 2010 November spectrum, 
but were not measurable for the 2011 December spectrum due to its 
lower S/N and resolution. 
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Fig. 5. Comparison of blue-end optical spectra obtained in 2010 November, 2011 July and 2012 July, after the resolution of the 
former has been degraded to match that of the latter two spectra (which themselves differ in resolution and wavelength coverage; 
Sect. 2). Prominent transitions are indicated as are the I-band magnitudes on the dates of the observations. Following the line 
identifications of Zickgraf et al. ( 1989 ), additional weak emission lines from low excitation metallic species such as Ca r, Cr n, Fe i, 
[S h], Si ii and V n are not shown for reasons of clarity. 



low excitation metallic lineal identified by previous studies are 
present in all four spectra, we find little indication of variability 
in their strength (cf. comparison of the 2010 November and 201 1 
July spectra; Fig. 5). Likewise, they demonstrate no appreciable 
evolution in width (FWHM~ 45 + 5kms _1 ; measured from the 
2010 November spectrum) since the observations of Zickgraf et 
al. ([1989) and imply much slower outflow velocities than the 
Balmer and He i P Cygni profiles do. Unfortunately, comparison 
of the line widths of the low excitation species between epochs 
is not instructive due to the lower resolution of our remaining 
spectra, but we see no appreciable differences in radial velocity 
in our data. 

We therefore arrive at a simple picture where the behaviour 
of high- and low excitation lines appears largely de-coupled, 
with the former showing dramatic variations in line profile and 
strength over timescales as short as months, while the latter re- 



3 Ca i, Cr ii, Fe i, Fe n, [Fe n] , [S n] , Si n, Ti n and V n. 



main largely unaffected. Moreover, this pattern appears qualita- 
tively consistent with observations made over the past half cen- 
tury (Sect. 1.1 and Table A.l). 

The observations available in the 6000-7200A and 8400- 
8900A windows likewise indicate that the overal spectral mor- 
phologies remain largely unchanged from previous epochs (Sect. 
1.1; Figs. 7 & 8). Specifically, and as expected, Ha and the 
Paschen series and a number of He i lines are strongly in emis- 
sion, although the I-band spectrum is dominated by emission 
from Ca n and O i. Other low excitation species present and in 
emission are Fe n, Si n and N i. Forbidden [O i] 6300A line is 
also in emission; while consistent with a [N n] identification we 
suspect the inflection in the blue flank of He is instead due to the 
run of absorption in the stellar wind (also seen in other sgB[e] 
stars such as R126; Zickgraf et al. 1986 1. Finally we identify TiO 
bandhead emission, as has previously been observed (Zickgraf et 
al. [19891 . 
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Unfortunately, the low excitation emission lines (e.g. [Oi], 
[On], Can and [Can]) employed by Rraus et al. (2010]) and 
Aret et al. (120121 1 to determine disc kinematics are either ab- 
sent, outside our coverage, or observed at too low a resolution 
to impart information on the velocity structure of the emitting 
region. However, in the absence of the shell-type emission line 
profiles that arise from obscuration of the stellar disc by an in- 
tervening circumstellar disc (cf. S65; Kraus et al. 2010 1 we infer 
the system is not observed edge-on. 

3.3. Photometric/spectroscopic comparison 

We may utilise our well sampled optical lightcurve to place the 
spectroscopic observations into context (Figs. 1 & 3); noting that 
unfortunately the latter sample a restricted photometric range 
(I~ 13.10 - 12.51 mag), albeit including local maxima and min- 
imcQ. Nevertheless, the data suggest that the changes in spectral 
morphology - and in particular the appearance of He n 4686A 
emission - is not uniquely correlated with photometric magni- 
tude (Table A. 1). This apparent lack of correlation is exemplified 
by the spectra obtained in 2010 November and 201 1 July, where 
the strength of the high excitation emission features in the spec- 
tra differ greatly (Figs. 5 and 6) but for which the I-band magni- 
tudes are identical to within ~0.07 mag. At the times of the two 
observations the extended photometric behaviour of the star was 
rather different, respectively exhibiting rapid low amplitude vari- 
ability and at the peak of a rapid (~ 60 d) ~0.35mag 'flare' (Fig. 
1), which might be supposed to be a cause of the different spec- 
tral morphologies. However, He n 4686A emission is present in 
both the 2011 and 2012 July spectra - obtained in local pho- 
tometric minima and maxima - implying a lack of a clear corre- 
spondence between lightcurve behaviour and spectral morpholo- 
gies. Nevertheless, our limited dataset is consistent with the pos- 
sibility that while He n 4686A emission may occur in any photo- 
metric state, a lack of emission is only associated with brighter 
states. Unfortunately, with only 5 quasi-contemporaneous obser- 
vations (Fig. 3) we are unable to draw any conclusions regarding 
correlations to the (V-I) colour; again, more frequent spectro- 
scopic sampling will be required to explore any such behaviour. 



4. The circumstellar environment of S18 

So how is S18 to be understood in the context of other sgB[e] 
stars? A defining property of such stars is the presence of a cool 
equatorial disc or torus, comprising neutral and molecular gas 
and warm dust. The presence of e.g. Ca n emission (Fig. 8) im- 
plies cool, dense conditions (~ 5000K, ~ 10 9 cm 3 ; Hamann 
& Simon 1988) which are shielded from the ionising flux of 
the central star. Likwise O i 8446A requires a dense transitional 
zone between H i and H ii regions (10 9 - 10 10 cm~ 3 ; Grandi 1980 ) 
for Ly/3 pumping to drive it strongly into emission (cf. Wdl-9; 
Clark at al. 120 13 al l and the discovery of Raman-scattering lines 
also requires neutral hydrogen to be present (Torres et al. I2012I) . 
The presence of TiO and CO bandhead emission (Sect. 1 . 1 and 
Fig. 7) likewise points to the presence of cool, dense molecular 
gas (Morris et al. [1996 Liermann et al. 120 10t . At longer wave- 
lengths the mid-IR properties of S18 suggest the presence of 
warm (~800K) dust with a probable crystaline silicate composi- 



4 van Genderen & Sterken (20021 report that the spectrum presented 
by Zickgraf et al. 0989} - in which Hen 4686A was absent - was 
obtained during a photometric maximum with V~13.3; comparable in 
magnitude to the peak in our lightcurve around JD~ 2455000. 
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Fig. 6. As Fig. 5 but spectra scaled to compare the stronger emis- 
sion features. 
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Fig. 7. As Fig. 5 but with the low resolution spectrum from 201 1 
December 8 - the only epoch for which coverage extends to Ha 
and the He 1 6678 and 7065 A transitions. Note the peak intensity 
of Ha is ~ 35 x continuum at this time. 



tion; both properties consistent with other sgB[e] stars (Kastner 
et al.EOOUIlSTOl Bonanos et al.|20T0]l. 

In these regards S 1 8 appears to conform to the stereotype of a 
sgB[e] star (Zickgraf et al. |19851 l, where the low excitation spec- 
trum is formed within - or in the case of the Lya and £1 pumped 
lines such as O i 8446A and Fe n 8490A - on the surface layer 
of a cool dense torus surrounding the central star. The optic al- 
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height and Qk(t) is the Keplerian angular frequency at radius r. 
Adopting M t ~ 19M and 7?, ~ 33R Q (e.g. Aret et al. l20l2l ) for 
r = R t we have fi^ ~ 0.8 d and H/r ~ c s /(r x Q K ) ~ 0.02 
and £2^ ~ 25 d and H/r ~ 0.06 for r = 10/?* (where c, is the 
sound speed). From these values we find f v i s ~ 20/(a/0.3) yr 
and ~ 62/(ff/0.3) yr respectively (where a is the viscosity pa- 
rameter and t V i S is proportional to the inverse of the disc temper- 
ature). Therefore the combination of the relatively low tempera- 
ture (in contrast to classical Be star discs) and large inner radius 
adopted results in a long viscous timescale and hence timescale 
of variability in response to changes in the bulk properties of the 
disc/torus. 

Despite this, both Nota et al. ( 119961 ) and Torres et al. (120121 1 
report line profile variability in the Fe n lines - one might expect 
rapid changes in the disc properties in response to changes in 
the ionising flux or, by analogy to classical Be stars, due to the 
impulsive injection of matter (Rivinius et al. 1998). 



Fig. 8. I-band spectrum of S18 from 2011 July 20 with major 
transitions indicated. 

IR continuum excess demonstrated by S18 (Zickgraf et al. 1989) 
would then result from a combination of thermal bremsstrahlung 
from the stellar wind and/or disc and thermal emission from the 
warm dust which also resides in the torus (e.g. Porter 2003 1. 

Recent spectroscopic observations suggest that the neutral 
gaseous component in sgB[e] stars is detached from the central 
star (e.g. Kraus et al. 120101 Liermann et al. 120101 & Aret et al. 
120121 1. Kastner et al. (2006 2010) arrived at a similar conclu- 
sion from analysis of the mid-IR properties of sgB[e] stars and, 
by analogy to lower mass AGB stars, postulated that the dust 
resides in long lived quasi-Keplerian disc surrounding a binary 
system which formed from mass lost in a previous common- 
envelope/mass transfer phase. Support for this geometry has 
been provided by interferometric observations of the galactic 
B[ e] star s CI Cam and MWC300 (Thureau et al. |2009l Wang e t 
al. l20T2l ) and the sgA[e] star HD 62623 (Meilland et al. [20T0b . 
which imply ring-like geometries for the circumstellar material. 

While our spectroscopic observations lack the wavelength 
coverage and resolution to address the kinematics of the 
disc/torus, the spectra of S18 presented by Aret et al. (2012) re- 
veal that the higher excitation [Ca n] lines are broader than the 
lower excitation [O i] lines, as anticipated for Keplerian rotation. 

Consequently, it is possible to see why the low-excitation 
component of the spectrum of S18 is, in terms of variability, es- 
sentially decoupled from the high-excitation component, ;/ the 
former arises in a largely static, quasi-Keplerian, viscous, de- 
tached disc while the latter is associated with an intrinsically 
variable star. Following the analysis of classical Be star discs 
by Okazaki et al. (119911) we may obtain an estimate of the vis- 
cous timescale of such a disc. We make the simplifying as- 
sumption that the gaseous component of the quasi-Keplerian 
disc is isothermal and adopt a representative disc temperature 
of Ta ~ 6OOOK0 and adopt representative inner disc radii of 
Rd ~ R* (i e. it extends to the stellar surface at R t ) and/?a ~ 10R, 
(motivated by the apparent truncation of the disc in LHA 1 15-S 
65; Kraus et al. [20TUl Aret et al. [20T2t . 

Under these conditions the viscous timescale in the disk, f v j s , 
is given by (a x (H/r) 2 x £2g-( r )) _1 , where H is the disk scale- 

5 The mid-point Kraus et al. (2010) derive from the strength of [Oi] 
emission lines, noting that as with classical Be stars the assumption of 
an iothermal disc is likely violated in practice. 



5. The nature of the central star(s) 

What is the nature of the hot, luminous central star in SI 8? 
No uncontaminated photospheric features are present in any 
published spectra so we are forced to infer stellar properties 
from the highly variable high-excitation emission-line spectrum. 
Following the suggestions of Zickgraf et al. ( 1 19891 ) and Nota et 
al. ( 119961 ), in Fig. 9 we plot examplar spectra from November 
2010 (henceforth the 'cool' spectrum/phase) and July 2012 (the 
'hot' spectrum/phase) against representative spectra of the early 
blue hypergiant (BHG) ( l Sco, the BHG/LBV P Cygni and 
a sample WN9h star, J013327.26+303909.1. As described by 
these authors, in both cool and hot states S18 shows gross 
similarities to the latter two objects. Specifically, (i) a simple 
spectrum comprising narrow Balmer and He i emission lines 
in the cool phase, with the former showing P Cygni profiles 
and (ii) the disappearance of P Cygni line profiles and the de- 
velopment/strengthening of He i, He ii and Nm emission in the 
hot phase; all characteristic features of the spectra of P Cygni- 
type LBVs and WN9h stars, respectively (e.g. Najarro 2001. 
Crowther & Smith 1997 ). Neither the hot nor the cool spectrum 
resembles that of the early BHG Sco and, by extension, field 
blue supergiants, because such stars do not support winds of suf- 
ficiently high mass loss rate to drive the Balmer series into the 
intense emission present in S18 (Table 1 and Clark et al. 12012b! 
and refs. therein). 

Under the assumption that the high-excitation emission line 
spectrum arises in a stellar wind, the blue edge of the P Cygni 
profiles in the Balmer series provides a lower limit to the wind 
velocity in the cool phase (~ 500 - 750kms~ 1 ). Such a value is 
broadly consistent with expectations for early-mid B supergiants 
(BSG)/BHGs, but much faster than expected for LBVs such as 
P Cygni (Table 1 and Clark et al. 2012b]). Conversely the mass 
loss rate (~ 3 x lO~ 5 M yr _1 ; Zickgraf et al. 1 19891 derived from 
line profile modeling of the Balmer series) is much higher than 
expected for BSGs and BHGs but consistent with known LBVs 
(Clark et al. [2012b The disappearance of P Cygni profiles in 
the WN9h spectrum is explicable in terms of the full ionisation 
of the stellar wind due to an increase in stellar temperature, al- 
though one might also expect such behaviour to be accompanied 
by changes in the stellar wind; note, for example, the greater 



6 Nevertheless, this pattern of wind properties is observed in the un- 
usual BHG Wdl-5; a star which is thought to be the result of binary 
driven mass loss and that is also thought to approach the Eddington 
limit (Clark etal. 120 13bt . 
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Fig. 9. Comparison of 2010 ('cool' state) and 2012 ('hot' state) spectra versus the Bl Ia + hypergiant Sco, P Cygni and the 
M33 WN9h star J013327. 26+303909.1 with major transitions indicated. Note that in the absence of a cool, toroidal circumstellar 
envelope, we would not expect the three comparison spectra to demonstrate the rich low-excitation metallic emission spectrum of 
S18. 

Table 1. Stellar parameters of comparison stars in Fig. 9 



Star 


log(L/L Q ) 


7eff 


M 




BC 


Ref. 






(kK) 


(lO^Moyr" 1 ) 


(km s" 1 ) 


(mag) 




C Sco (Bl Ia + ) 


5.93 


17.2 


1.55 


390 


-1.6 


1 


P Cygni (Bl Ia + /LBV) 


5.85 


19.0 


30.0 


185 


-1.8 


2 


Mean WN9h 


5.7 


27.6 


18.5 


450 


-2.9 


3 



Tabulation of the bolometric luminosity (L), effective temperature (r e ff X mass loss rate (M), terminal wind velocity (v m ) and bolometric correction 
(BC) for the comparison spectra in Fig. 9. References are 'Clark et al. d2012bt . 2 Najarro ([200 1 > and 3 Crowther & Smith ( 1997 (; the value for the 
final entry is the mean of the modeling results for the four WN9h stars analysed by the latter authors. The broadly comparable luminosities reveal 
that the changes in spectral morphology are primarily the result of changes in temperature and wind parameters. 



wind velocity of WN9h stars with respect to P Cygni-type LBVs 
(for similar mass loss rates; Table 1). Under such a picture the 
spectrum obtained in 201 1 December, which demonstrates weak 
He n emission and P Cygni profiles in the higher Balmer transi- 
tions, would represent an intermediate state. 



An obvious interpretation of these spectroscopic data is that 
S18 transitions directly between P Cygni-type LBV and WN9h 
spectral types. However there are several caveats to such an hy- 
pothesis. Firstly, Zickgraf et al. (1989) presented a simple model 
atmosphere fit to the UV-IR spectral energy distribution of S18 
and, assuming a temperature of ~25000K, a bolometric correc- 
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tion of -2.5mag and a distance to the SMC of 60kpc, inferred 
a luminosity of log(L/L ) ~ 5.6 in the cool phase. The lower 
bolometric correction derived for P Cygni suggests a downwards 
revision to log(L/L Q ) ~ 5.3; significantly lower than in the com- 
parison stars listed in Table 1 (although low luminosity examples 
are known - e.g. Clark et al. 2009a i. 

Secondly, in both cool and hot phases the emission line spec- 
trum is unprecedentedly strong with respect to P Cygni-type 
LBVs and WN9h stars; this being particularly noteworthy for 
Hen 4686A (Fig. 9: Nota et al. [19961 Moreover, Hen 4686A 
was seen to be double peaked at one epoch, with the blue peak 
of greater intensity (with corresponding asymmetric profiles in 
the Balmer series and Hei; Nota et al. 11 9961 1 . Such asymme- 
tries are observed in other B[e] stars (e.g. Zickgraf 2003 1 and are 
also common in classical Be stars (Okazaki ll991l l; in all cases 
heralding a departure from spherical symmetry for the circum- 
stellar environment, as we infer for SI 8. This is in contrast to 
the spherical geometry typically assumed for the the winds of P 
Cygni-type LBVs and WN9h stars (although Nota et al. ( 119961 
report similar asymmetries for a subset of the Ofpe/WN9 stars 
they studied). Another unexpected spectral peculiarity in S18 is 
the absence of P Cygni profiles in the He i emission lines in the 
cool spectrum when clearly present in the Balmer series. 

Finally, the X-ray detection (Lx ~ 3 x 10 33 erg s _1 ; Antoniou 
et al. 120091 Bartlett et al. in prep.) strongly suggests a binary 
nature for SI 8. The canonical X-ray luminosity for a single 
massive star (Lx ~ 10~ 7 Lb o i; Berghoefer et al. 119971 1 implies 
Lx ~ 1.6 x 10 32 erg s _1 for SI 8; over an order of magnitude 
lower than observed (Sect. 1.1) and hence requiring an additional 
source of X-rays to be present within the system. Therefore in 
addition to a contribution from the circumstellar disc/torus one 
might also suppose additional spectroscopic contamination from 
a binary companion; we return to this below. 

5.1. S18asaLBV 

With these caveats in mind, the simplest explanation for the 
spectrosopic and photometric properties of S18 remains that it 
hosts an LBV. Indeed, a physical link between sgB[e] stars and 
LBVs has been proposed on the basis of photometric variability 
in a number of starsj. 

S18 appears to reside in a region of the H-R diagram prone 
to pulsational instabilities (e.g. Cantiello et al. I2009I Saio l201 11 1 
and clearly shows the presence of the rapid microvariability 
present in many LBVs and luminous super-/hypergiants (the a- 
Cygni variables; van Leeuwen et al. [1998 □ but the rapid, high- 
amplitude variability it demonstrates appears at variance with 
typical canonical LBV excursions (e.g. Spoon et al. [T9941PI . 

7 T he binary R MC 4 (Zickgraf et al. [1996V LHA 120-S 22 (Shore 
et al. 11990119921 van Gendere n & Sterken [Tg9"9l Sterken I2UTT1 
LHA 120-S 134 (Stahl et al. [T9841 van Genderen vanGOl) and 
J013459.47+303424.7 (Clark et al. I2012at are all reported as photo- 
metrically variable at UV-optical wavelengths, albeit from considerably 
sparser datasets than considered here. 

8 Following the analysis of the blue supergiant HD 50064 by Aerts 
et al. (20101, an identification of these variations with strange mode os- 
cillations would appear an interesting hypothesis, especially given the 
long-term changes in both amplitude and timescale of this behaviour. 
However, the likely progenitor mass of S18 (~ 20 - 25M B via compar- 
ison to the isochrones of Groh et al. 12013) is lower than that expected 
for stars encountering this instability (~ 40M o ; Aerts et al. 20101. 

9 These consist of (i) microvariations, occuring over weeks to months 
and of low amplitude (< 0.2 mag), (ii) LBV excursion, occuring over 
years to decades and typically of ~ 0.5 - 2.0 mag amplitude and which 



Nevertheless such rapid variability is not completely unknown, 
being present in a handful of local LBVs and also the ex- 
tragalactic SNe imposters, presaging both giant eruptions (rj 
Cannae, HD5980 and S N 1954J; Frew I20U41 Smi th & Frew 
1201 11 Breysacher 1997 and Tammann & Sandage 1968) and 
bona fide SNe (SN 2010mc and possibly SN2009ip; Ofek et al. 
120131 and Sect. 6), although in the case of SN2000ch a decade 
of such behaviour has yet to lead to either outcome (Pastorello 
et al. l20T0l l. S18 may be identified with SN2000ch in this re- 
gard, although the amplitude of variability is significantly lower, 
most closely resembling SN 1954J in the 4 years leading up to 
its eponymous outburst (Tammann & Sandage 1968; see also 
SN2009ip (Sect.6)). 

The colour behaviour of S18 also differs from that of LBV 
excursions; leading to the system becoming bluer in the ris- 
ing branches of the lightcurve and showing much greater vari- 
ance in comparison to canonical examples (cf. Fig. B.l, de- 
rived from data presented by Spoon et al. 119941 . The colour 
behaviour of normal LBVs - which become redder as they 
brighten - arises from the simultaneous expansion (contraction) 
and cooling (warming) of the photosphere as the star transits 
from WN10-1 lh to cool supergiant in 'outburst' and back again 
(e.g. AG Car; Groh et al. |20U91 Groh priv. comm. 2012) and 
clearly cannot explain the behaviour observed here. Indeed the 
spectral evolution of S18 occurs over a temperature range for 
which the change in (V-I) colour index would be minimal (cf. 
Clark et al. l2005al l. 

We highlight that not all LBVs show this behaviour - for in- 
stance the SNe imposter SN 2000ch shows no colour/magnitude 
correlations during its rapid high amplitude flares (Fig. B.l; 
Pastorello et al. 1201 Oi l. An additional complication with inter- 
preting the behaviour of S 18 is that the presence of a substan- 
tial circumstellar envelope - and a probable binary campanion 
- may also contribute to the photometric variability; indeed the 
modeling of Zickgraf et al. ( 1989) suggests a continuum excess 
extending to optical wavelengths. 

We may advance two different scenarios to explain the pho- 
tometic variability, noting that they are not mutually exclusive. 
Firstly, the episodes in which S 1 8 rapidly brightens and becomes 
bluer may be due to an increased contribution from the hot 
central star relative to the cool circumstellar disc/torus (and/or 
stellar wind) in the combined emergent spectrum. The lack of 
colour changes associated with 'quiescent' and fading branches 
of the lightcurve is more difficult to understand, but may rep- 
resent quasi-simultaneous variability in both emission compo- 
nents. Moreover, we are unable to advance a physical mecha- 
nism for such stellar 'pulsations', although this is true for LBVs 
in general; for instance the pulsations Aerts et al. (20101 attribute 
to strange mode oscillations in HD 50064 are of substantially 
lower magnitude (~ 0.2mag) to those observed here. 

Alternatively the variability may be due to changeable cir- 
cumstellar extinction, with comparable behaviour seen in both 
classical Be (de Wit et al. 120061) and Herbig Be stars (Lamers et 
al 1 1999I >. In classical Be stars this is observed for edge-on sys- 
tems, where the gaseous disc projected against the star reduces 
continuum emission at short wavelengths, but contributes con- 
tinuum emission at predominantly longer wavelengths from the 
portion of the disc not projected against the star, such that the 



are typically associated with colour changes in the sense that the star 
becomes redder as it brightens and (iii) giant eruptions, in which the star 
brightens by > 2.0 mag over timescales ranging from days to years and 
in which the bolometric luminosity is not conserved (e.g. Humphreys 
& Davidson [T994l Lamers et al. lI99"8l . 
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system becomes fainter and redder as the disc optical depth in- 
creases. A variant of this scenario is assumed for Herbig Be stars, 
where variable stellar extinction is attributable to a clumpy, dusty 
circumstellar envelope. 

However, a significant problem with such an hypothesis 
is the location of the obscuring material and the rapidity of 
changes. A comparison of the line profiles of S 18 to e.g. those 
of the sgB[e] star LHA 1 15-S 65 suggest it is not seen edge-on 
(Zickgraf et al 1989, Rraus et al. 2010); therefore it appears un- 
likely that the circumstellar disc/torus could be responsible as is 
assumed for both classical and Herbig Be stars. Likwise, while 
the stellar winds of LBVs are known to be clumpy (Davies et al. 
2005 ), such clumps are not expected to be sufficiently dense to 
provide the requisite obscuration of the stellar disc. 

Irrespective of the physical mechanism, a key observational 
finding of this research is that there is no clear correlation be- 
tween photometric and spectroscopic variability; for example the 
hot and cool phase spectral states in 2010 November and 2011 
July differed by only ~ 0.07mag. Naively applying the respec- 
tive bolometric corrections implied by the spectral morphologies 
(Table 1 ) would imply that the luminosity of S 1 8 had increased 
by > OAAdex over this period; however such a conclusion is 
tempered by the potentially highly variable continuum contribu- 
tion from the stellar wind/circumstellar envelope at these wave- 
lengths. Bearing this important caveat in mind, changes in bolo- 
metric luminoity are not unprecendented in bona fide LBVs, 
having been observed in e.g. AFGL 2298 and Romano's star 
(Clark et al. |2009bl Maryeva & Abolmasov l2()T2l . Moreover, 
while S18 has not been observed to transition to the cool hy- 
pergiant phase that characterises many LBV excursions, the 
range of spectroscopic variability it exhibits (WN9hoP Cygni- 
type LBV) mirrors the behaviour of Romano's Star (Maryeva & 
Abolmasov l20l2l . 

Finally, we note that the mid-IR properties of S18 appear 
to exclude the presence of the cool dust that characterises the 
circumstellar ejecta of many LBVs (Clark et al. 2003. Bonanos 
et al. l20T0b . 

5.2. S18as a binary system 

The presence of highly variable He n 4686A emission in the 
spectrum of S18 has long been attributed to binarity (Shore et al. 
flWTl Zickgraf et al. H98^ . Furthermore Naze et al. (l20T2l find 
that in the absence of a binary companion LBVs are intrinsically 
X-ray faint, strongly suggesting that our X-ray detection of S18 
indicates the presence of an unseen binary companion (Sects 1 . 1 
and 4.0). Two options therefore present themselves; accretion 
onto a compact object or wind collision leading to shocked ma- 
terial in a massive binary. 

Several examples of X-ray binaries with sgB[e]-like mass 
donors are known. SS433 has a complex, dusty circumstellar en- 
vironment (Clark et al. 2007 1 and demonstrates highly variable 
He n 4686 A emission (Fabrika & Bychkova 1990), but is several 
orders of magnitude brighter at X-rays than S18 (Cherepaschuk 
et al. 120051 CI Cam and IGR J16318-4848 appeal- to provide 
better comparisons; while distances to both stars are uncertain, 
at the lower end of the ranges quoted in the literature their qui- 
escent fluxes are comparable to that of S18 (Bartlett et al. 120 121 
Barraganet al. 120091 . 

Spectroscopically, CI Cam also bears a close resemblance, 
even to the extent of He n 4686A emission being present in qui- 
escence (e.g. Clark et al. 1999 Hynes et al. 2002). Nevertheless, 
the long term photometric and spectroscopic behaviours of CI 



Cam and S18 are very different, with the former much less vari- 
able and demonstrating consistently weaker He n 4686A emis- 
sion (Clark et al. l2000l Hynes et al. 120021 . Indeed, He n 4686A 
only approached the strength observed in S18 when CI Cam ex- 
hibited an X-ray flare in 1998 April. The outburst was replicated 
at optical wavelengths, where CI Cam brightened by over 2.5 
magnitudes (Clark et al. 2000); to date no such behaviour has 
been observed for SI 8. 

Conversely, there are a number of observational precedents 
that suggest that S18 could comprise a massive interacting bi- 
nary system. Firstly, both GG Carina and VFTS698 are spec- 
troscopically confirmed sgB[e] binaries with apparent circumbi- 
nary discs/toroids (Gosset et al. 119851 Rraus et al. 120131 Dunstall 
et al. 120121 ). while Kastner et al. (1200611201 01 suggest such a con- 
figuration from their modeling of the mid-IR excess of sgB[e] 
stars. Secondly, the X-ray luminosity of S 18 is well within the 
range observed for colliding wind systems (~ 10 32 - 10 35 ergs~'; 
Gagne et al. 120121 ) and is comparable to that of the sgB[e] star 
Wdl-9; a massive compact interacting binary within the star- 
burst cluster Wdl (Clark et al. 2008. 2013a). 

Finally, variable He n 4686A emission is a key observational 
signature of a number of massive binary systems, where it is 
the result of an additional contribution to the ionising flux from 
shocks arising in a wind-wind collision zone. Under such a 
scenario we might explain the increase in strength of the He n 
4686A line in tandem with the apparent transition to a hot, 
WN9h phase as a result of the increase in wind velocity expected 
for a star of such a spectral type (Table 1), which in turn would 
lead to a stronger wind collision shock. 

In this regard three LBVs are of interest, rj Carina is an 
eccentric, long period (P a> ~ 5.52yr), colliding wind binary 
that shows orbital modulation in the He n 4686A line (Steiner 
& Daminelli 2004; Teodoro et al. 120121 ) - albeit much weaker 
that than seen in S18 - and strong X-ray emission (maximum 
L x ~ 1.3xl0 35 erg s _1 ; Ishibashi et al. [19991 

The second, HD 5980, is a close, eccentric, eclipsing bi- 
nary (f or b ~ 19.26d; Sterken & Breysacher 11997) compris- 
ing an LBV primar y (WN 3oWNl lh/ B1.5 I a + ; Niemela IT9881 
Koenigsbergeret al. 1994, Drissen et al. 2001 Heyardi-Melayari 
et al. 1 19971 ) and a WN4 secondary (Breysacher et al. 1982 ) with 
a possible third 04-6 component in a wide orbit (Foellmi et al. 
2008). Highly variable Hen 4686A emission has been observed 
both before and after its 1993-5 outbursts (e.g. Moffat et al. |1998l 
Breysacher & Francois 2000. Koenigsberger et al. [2010 >. This 
is likely the result of wind/wind collisional effects, although an 
intrinsic contribution from one or both stellar winds cannot be 
excluded. Unlike SI 8, this emission does not appear to be ab- 
sent at any (orbital) phase. X-ray emission is present and of a 
comparable magnitude to S18 (Lx ~ 1.7xl0 34 erg s _1 ; Naze et 
al. l2TJ02l 120071 . 

The final star is the SNe imposter SN 2000ch, which also 
demonstrates an asymmteric doubled lined He n 4686A line in 
quiescence (Pastorello et al. 2010). These authors draw direct 
comparison between SN 2000ch, S18 and indeed HD 5980 on 
this basis, suggesting that all three could be binaries and we note 
that both SN 2000ch and HD 5980 also exhibit the rapid photo- 
metric variability seen in SI 8, albeit of much greater amplitude. 

Therefore the multiwavelength properties of S18 are at least 
consistent with other luminous colliding wind binaries contain- 
ing LBVs. Given the unprecented strength and variability in the 
He ii 4686A line it would be of considerable interest to search 
for periodic variability indicative of orbital modulation, espe- 
cially since the inclination of S18 (Zickgraf et al. 1989) dis- 
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favours identification of binary reflex motion and no evidence 
of eclipses or elipsoidal modulation is present in the lightcurve. 
Moreover, higher S/N X-ray observations in order to distinguish 
between relativistic-accretor and colliding-wind scenarios would 
be of considerable interest. 

6. Concluding remarks 

Photometric and spectroscopic observations of S18 over the past 
half century demonstrate that it is a highly dynamic and complex 
system; at odds with the historic picture of sgB[e] stars as essen- 
tially static objects. The cool circumstellar disc/torus shows little 
spectroscopic variability over the course of this period; appear- 
ing to be essentially decoupled from the central stellar source, 
which is delineated by the high excitation emission component 
of the composite spectrum and is highly variable. 

Our observations reveal that, trivially, S18 satisfies the 
eponymous classification criteria of LBVs. However, this does 
not imply that the physical causes of its variability are the same 
as those of canonical S-Doradus type objects such as AG Car. 
Specifically, the (assumed) 'stellar' spectrum - and the speed 
of its evolution - is more extreme than canonical examples. 
Moreover, the photometric variability of S 18 is difficult to inter- 
pret in the framework of known LBVs, although the rapidity of 
the 'flaring' does appear to mirror that of some 'SNe imposters'. 
The colour behaviour is likewise atypical and, surprisingly, we 
see no apparent correlation between spectroscopic and photo- 
metric behaviours, raising the possibility that the bolometric lu- 
minosity of S18 varies over time. 

The detection of X-ray emission adds to the weight of evi- 
dence that S 18 hosts an unseen binary companion, with observa- 
tions apparently favouring a colliding wind system over a rela- 
tivistic accretor, although the latter cannot currently be excluded. 
In this respect the combination of a powerful wind at relatively 
modest luminosity inferred by Zickgraf et al. (1989 1 is of inter- 
est, since it mirrors the properties of the BHG Wdl-5, which is 
thought to have evolved through a phase of binary-driven mass- 
loss (Clark et al. I2013bl >. Model atmosphere analysis of S18 to 
determine if it shares the anomalous chemistry of Wdl-5 would 
be welcome, but is currently premature given its aspherical cir- 
cumstellar envelope and the possibility of spectral contamination 
resulting from binarity. 

These observational properties highlight the importance of 
S18 for a number of reasons. Firstly it strengthens the physical 
link posited between sgB[e] stars and LBVs and also the possi- 
bility that the B[e] phenomenon is linked with binarity; specifi- 
cally that the circumstellar/binary disc is the direct result of bi- 
nary evolution (e.g. Zickgraf 2003 Kastner et al. 120101 Kraus 
et al 1201 3b . The presence of a detached, quasi-Keplerian vis- 
cous disc/torus is attractive since the decadal viscous timescale 
of such a structure would naturally explain the relative lack of 
variability in its spectral signatures in comparison to the cen- 
tral star(s). Indeed, to the best of our knowledge, and in con- 
trast to classical Be stars, disc-loss and/or reformation has not 
been observed in sgB[e] stars. This would be explicable if the 
viscous (and hence dissipation/formation) timescales of discs 
around sgB[e] stars is significantly longer than that of classical 
Be stars due to their greater physical extent and lower tempera- 
tures. 

Recent observational studies suggest both a high binary frac- 
tion amongst massive stars and that the majority of these systems 
will in turn interact (Sana et al. 2012 1. Given this, the sugges- 
tion that sgB[e] stars may be currently undergoing, or have re- 
cently undertaken, significant binary-driven mass-loss/transfer is 



important, since it would allow for their identification with this 
brief evolutionary phase. Other examples of this phenomenon 
would then be the sgB[e] star Wdl-9 (Clark et al. 12013a) and 
RY Scuti (e.g. Smith et al. 1201 11 1; the latter system is of partic- 
ular relevance since the authors speculate that the mass transfer 
in this system is mediated by episodic, LBV-like events. 

Finally, the presence of a massive detached disc/torus asso- 
ciated with an apparent LBV raises issues regarding the nature 
of its demise, given recent observational and theoretical pre- 
dictions that such stars can serve as the immediate progenitors 
of core-collapse SNe (e.g. Gal-Yam & Leonard 2009 Groh et 
al. 120131 1. Indeed the latter authors find this occurs for stars in 
the range log(L bo i/£©) ~ 5.3 - 5.6 (M initia i ~ 20 - 25M ); di- 
rectly spaning the luminosity range we infer for S18 (Sects. 5 
and 5.1). Following from this it is immediately tempting to draw 
parallels between S18 and the B3 la progenitor of SN1987A 
(Walborn et al. 119891 1, which was of similarly modest luminos- 
ity (log(Lt,oi/£o) ~ 5.1 ±0.1), must have been associated with 
a circumstellar disc and for which a binary interaction/merger 
scenario has been proposed that would have left the progenitor 
resembling a sgB[e] star in its immediate aftermath (e.g. Morris 
& Podsiadlowski 12509b . 

Another intriguing comparator is SN2009ip. Initially clas- 
sified as a SNe-imposter, between 2009-11 it demonstrated 
rapid photometric variability of a similar, but more extreme na- 
ture than S18 ( Sect . 5.1; Smith et al. I20TU1 Foley et al. I20TT1 
Pastorello et al. 120131 1. During this period, its spectrum was also 
not dissimilar to that of S 18 in gross terms, being dominated by 
Balmer series He i and Fe n line emission in the blue and Paschen 
series, Ca n and O i emission in the red, although outflow veloci- 
ties were substantially greater than implied for S 1 8 and a greater 
mass (~ 50 - 80M o ) also inferred. In 2012 it displayed two erup- 
tions; the nature of the first event is uncertain, with core-collapse 
SNe, pair-production instability and binary merger all being pro- 
posed, while the second is assumed to result from the interaction 
of rapidly expanding ejecta with existing circumstellar material 
(Fraser et a l. 12013] Mauerhan et al.l20T3] Pastorello et al. 120131 
Prieto et al. l20T2l S oker & Kashi l20T3l) . Both Smith et al. (120131 ) 
and Levesque et al. (2012) conclude that dust was present in the 
circumstellar material, with the latter authors suggesting a disc- 
like geometry and a possible origin in binary-driven mass-loss 
(although an origin in the recent photometric outbursts is also 
viable). As such, the possibility that S 1 8 represents a lower mass 
analogue of this highly unusual energetic transient is particularly 
exciting. 
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Table A.l. Optical spectroscopic summary 



Date 


I-band 


Comments 


Reference 


< 1955 




Ha, H/3, Hy and He n 4686A in emission 


Lindsay ( 1955) 


1967 July 14 




Balmer emission to He. He n 4686A absent 


Sanduleak d!978t 


1972-74 




Ha, H/3, Hy, He in emission. No mention of He n 4686A 


Azzooardi et al. (119751 


1977 September 17 


- 


Balmer emission to He. Hen 4686A emission comparable to H/3 


Sanduleak (1978) 


1978 January 27 


- 


H/3 to H13 P Cygni emission. Hen 4686A and Fen (multiplet 24) in emission 


Azzopardi et al. dl981) 


1978 July 16 




Ho- to HI 1 and He n 4686A in emission 


Shore etal. (1987) 


1981 July 13 


- 


Balmer series and He n 4686A in emission 


Shore & Sanduleak|l982| 


1982 March 22 


- 


Balmer series and He n 4686A in emission 
Hen 4686A -50% weaker than 1981 July 13 


Shore & Sanduleak[l982] 


1983 Oct 10 


- 


Ho- to HI 1 and He n 4686A in emission 


Shore et al. dl987t 


1983 December 17 


- 


Ha to HI 1 and He n 4686A in emission. 

P Cygni profile for higher Balmer lines but not Ha or He n 4686A 


Shore et al. (1987) 


1987 Nov. 7 




Ha to H(5 and He i 5876A in P Cygni emission. He n 4686A absent 


Zickgraf et al. (1989 1 


1991 September 17-20 




Ha to US and Hei in emission. Double peaked Hen 4686A 


Nota et al. (1996) 


2000 October 11-14 


12.52 


H/3 to HiS in emission. Hen 4686A absent 


Massey & Duffy (200U 


2000 October 13 


12.52 


Ha, H/3 and He i 6678A in P Cygni emission. He n 4686A absent 


Torres et al. (2012) 


2001 November 24 


12.66 


Ha, H/3 and Hei 6678 A in emission. Hen 4686A absent 


Torres et al. (2012) 


2005 December 10 


12.69 


Ha, H/3, He i 6678 A and Hen 4686A in emission 


Torres et al. (2012) 


2007 October 3-4 


13.10 


Ha, H/3, He i 6678 A and Hen 4686A in emission 


Torres et al. (2012) 


2008 November 13 


12.79 


Ha, H/3, He i 6678 A and Hen 4686A in emission 


Torres et al. (2012) 


2010 November 24 


12.55 


H/3 to He in P Cygni emission, He n 4686A absent 


This work 


2011 July 11 


12.62 


Ha to H5, He i 6678A and 7065 A and He n 4686A in emission 
He ii 4686A stronger than Hy 


This work 


2011 December 8 


12.51 


H/3 to He in emission, He n 4686A weakly in emission 
H/3 single peaked, Hy to He P Cygni profiles 


This work 


2012 July 07 


12.92 


H/3 to He in (single peaked) emission, He n 4686A stronger than Hy 


This work 



I-band magnitudes quoted are from OGLE-II,-III and -IV (Sect. 2.1). Note that emission lines are single peaked unless otherwise noted. Lindsay 
(1955), Sanduleak d!978t . Azzopardi et al. (1975), Shore & Sanduleak d!982t . Shore et al. d!987t and Torres et al. d2012) do not present all the 
spectra described in these works. Additionally Azzopardi & Breysacher (1975) detected He n 4686A in emission in 1977 October via an objective 
prism survey. 



Appendix A: Summary of historical and current 
data for S1 8 

Appendix B: Colour magnitude behaviour of LBVs 
and SNe imposters 
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